An entangled two photon source using biexciton emission of an asymmetric quantum 

dot in a cavity 
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A semiconductor based scheme has been proposed for generating entangled photon pairs from the 
radiative decay of an electrically-pumped biexciton in a quantum dot. Symmetric dots produce po- 
larisation entanglement, but experimentally-realised asymmetric dots produce photons entangled in 
both polarisation and frequency. In this work, we investigate the possibility of erasing the 'which- 
path' information contained in the frequencies of the photons produced by asymmetric quantum 
dots to recover polarisation-entangled photons. We consider a biexciton with non-degenerate inter- 
mediate excitonic states in a leaky optical cavity with pairs of degenerate cavity modes close to the 
non-degenerate exciton transition frequencies. An open quantum system approach is used to com- 
pute the polarisation entanglement of the two-photon state after it escapes from the cavity, measured 
by the visibility of two-photon interference fringes. We explicitly relate the two-photon visibility to 
the degree of Bell-inequality violation, deriving a threshold at which Bell-inequality violations will 
be observed. Our results show that an ideal cavity will produce maximally polarisation-entangled 
photon pairs, and even a non-ideal cavity will produce partially entangled photon pairs capable of 
violating a Bell-inequality. 

PACS numbers: 78.67.Hc 12.20.Ds 03.65.Yz, 42.50.Ct 42.50.Lc, 



I. INTRODUCTION 

Recent proposals for quantum communication |^| 
and quantum information protocols provide a signifi- 
cant incentive to develop practical single photon sources 
and entangled two photon sources. The first requirement 
for such sources is that the emission times of the photons 
be periodic with a precisely defined clocked frequency. 
Exciton recombination in electrically or optically excited 
quantum dots is a candidate system for such sources. 
In this paper we will discuss an entangled two photon 
source based on recent experiments in self assembled in- 
terface quantum dots j|, || . A proposal for producing 
entangled photon-pairs on demand based on biexciton 
emission from a quantum dot was recently presented by 
Benson et al. [Q] . 

A pair of excitons confined in a quantum dot form a 
bound state known as a biexciton. The decay of the biex- 
citon proceeds by consecutive single electron hole recom- 
bination processes. This is estabished experimentally by 
the temporal correlation of the biexciton emission and 
the exciton emission; time resolved photoluminescence 
measurements show the exciton photon to be emitted 
after the biexciton photon A similar time resolved 
study of the polarisation of the emitted photons shows 
that there are two decay paths, and it has been shown 
that they are coherent with one another J7j. While the 
biexciton photon and the exciton photon emitted in each 
decay path have the same linear polarisation, the polar- 



istion in different decay paths are orthogonal. If these 
decay paths were indistinguishable then this would be a 
good candidate for an entangled two photon source. Un- 
fortunately small asymmetries in the physical geomery of 
the dots makes the two paths distinguishable, since the 
asymmetry of the dot breaks the degeneracy of an inter- 
mediate exciton level enabling the two paths to be dis- 
tinguished by frequency. The effect of asymmetry on the 
spectrum of excitons in dots was observed experimentally 
in dots formed by monolayer fluctuations in a GaAs 2D 
quantum well || and has been addressed theoretically || . 
It has also been observed experimentally in CdSe/ZnSc 
dots lliol and in self assembled GaAs/InGaAs dots [O. 
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In Fig. 1(a) we indicate the possible decay paths from a 
single biexciton level through two non degenerate exci- 
ton levels to the ground state of the dot. The first decay 
path corresponds to the emission of a biexciton photon 
with linear polarisation in the x-direction at frequency 
u>i, followed by the emission of the exciton photon, with 
the same polarisation, at frequency UJ2- In the second 
decay path the biexciton emits a y-polarised photon at 
frequency W3 followed by the exciton emission, also with 
y-polarisation, at frequency llh. 

The state of the emitted photon pairs may then be 
written as 

|V>i) = (|x,wi;x,w 2 ) + |y,w 3 ;y,w 4 ))/\/2, (1) 

where the notation indicates the mode (polarisation 
and frequency) occupied by each photon of the pair, 
I photon 1; photon 2), with the order reflecting the order 
of emission. It has been established experimentally that 
the weights of the kets are equal [jll). In contrast, we 
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wish to produce a state of the form 

\ip 2 ) = |x,om;x, ub) + |y, ojA;y, vb)/V2, 



(2) 



and we call such a state polarisation entangled, since the 
entanglement is only in the polarisation degree of free- 
dom. This is in contrast to the state given by Eq. (|l|) 
which is entangled in both polarisation and frequency. 
The important difference between states \ipi) and \ip2) 
is that the second ket in ^2) may be rotated into the 
first ket using linear optical elements such as half-wave 
plates (HWP) and polarising beam splitters (PBS), and 
vice-versa, whereas this is not possible for the two kets 
written in state \ipi). Thus, for instance, Bell-inequality 
measurements and two-photon interference experiments 
may be performed with realtive ease using \tpi) but not 
1^2)) an d this translates to a technological setting in, for 
instance, quantum key distribution. 

The problem of producing frequency-and-polarisation 
entangled states akin to \ipi) has been considered for 
photon pairs produced by spontaneous parametric down- 
conversion in a nonlinear crystal |l2|. In this case pho- 
tons are also entangled both in polarisation and fre- 
quency, though the frequency entanglement is more com- 
plicated. The frequencies for the two emitted photons are 
constrained by energy conservation, so that their sum 
must equal the frequency of the absorbed pump pho- 
ton. Since this single constraint does notdetermine the 
frequencies of the two emitted photons uniquely, each 
photon of the pair may be emitted over a wide range 
of frequencies determined by the spectrum of the pump 
pulse and the phase matching requirement (which is an 
expression of momentum conservation). Thus the photon 
pair is entangled in its frequency degree of freedom. 

A resolution to this problem, presented and experi- 
mentally implemented in is to pass the signal and 
idler beams back through the crystal, but with the po- 
larisations rotated through n/2, with the result that the 
two ways in which the photons can be emitted with cor- 
related polarisation arc not distinguished by frequency. 
This scheme in |Q docs not directly translate to the 
case of biexcitonic emission, but we nevertheless wish to 
remove the spectral dependence from the entanglement 
in state l^i), so the aim of this paper is to present and 
analyse a proposal to accomplish this for the biexciton 
entangled photon source. 

In this paper we demonstrate that the frequency may 
be disentangled from the polarisation by placing the dot 
in an external cavity with suitably chosen cavity-exciton 
coupling strengths and cavity mode frequencies. We will 



show that the external cavity can erase the "which-path" 
information contained in the frequency components of 
state The external cavity is used to control both 

the spectral and spatial mode structure of the emitted 
photons to enable the entanglement to be demonstrated 
in an interferometer. A similar idea using waveguides for 
SPDC has been proposed by Banaszek et al 13 1. We note 
that the original proposal for the two-photon source || 
includes the external cavity, but its presence is only to in- 
crease the outcoupling efficiency, and only brief mention 
is made of its effect upon the spectral emission properties 
of the emitted photons. 

The next part of this paper begins by defining a Hamil- 
tonian for a four level system interacting with optical 
cavity modes. A master equation is developed in the 
third section to deal with photons leaking from the cav- 
ity and into some measurement apparatus, as well as to 
account for decoherence events such as photon loss. In 
the fourth section, we discuss some operational defini- 
tions to quantify the entanglement of the photons pro- 
duced, such as two-photon visibility and Bell-inequality 
violations, with the aid of which we judge the efficacy of 
the cavity in restoring the polarisation entanglement. We 
then provide some results in the fifth and sixth sections 
showing that an ideal cavity does establish maximally 
entangled photon pairs, and numerical results showing 
how sensitive the resulting state is to imperfections in 
the system parameters. We then provide some heuristic 
analytic results in the discussion which explain the nu- 
merical results, as well as comment on implications for 
experiments, and finally conclude the paper. 



II. SYSTEM HAMILTONIAN 



Figure 1(a) shows the energy levels and available dipole 



transitions for the biexciton-cavity system. The biexci- 
ton state is given by \XX), \Xx) and |Xy) are the inter- 
mediate excitonic states in the x and y polarisation decay 
paths respectively and |G) is the dot ground state. The 
cavity is assumed to support pairs of degenerate x- and y- 
polarised modes at frequencies u>a and lob- In our model, 
we do not include coupling between for instance the cav- 
ity mode |o>a,x) and the transition |G) <-> |Ax) which is 
valid assuming the detuning between them is much larger 
than the cavity-exciton coupling strength, which is the 
case for this system. The system Hamiltonian, H SYS , un- 
der the rotating wave and dipole approximations |^4[ [lf| 
is then 



H sys = uo\XX){XX\+u 2 \Xx)(Xx\+uj 4 \Xy)(Xy\+UA{n x ,u A +ny,u, A )+UB 

+ l -( qi \Xx)(XX\al WA +q 2 \G)(Xx\al WB + q 3 \Xy)(XX\a$ t „ A + ^G) (Ay|a^ _ H .c), (3) 
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FIG. 1: (a) Energy level diagram and available transitions for the quantum dot and cavity system, (b) Spectrum of exciton 
transitions (dotted line) and cavity modes (solid line) indicating the relevant frequencies for the interaction Hamiltonian. 



where a,j and hj = ajoj are respectively the photon anni- 
hilation operator and photon number operator for mode 
j, and for convenience we take H = 1. We transform to 
an interaction picture defined by Ho = ^N, where 

N = 2\XX)(XX\ + \Xx)(Xx\ + \Xy)(Xy\ 



is the number of excitations in the system. The inter- 
action Hamiltonian, H — e lHat H sys e~ lHot — Hq, is given 
by 



H = - £\Xx)(Xx\ - (£ + A)|Xy)(Xy| + (£ + 6 A )(n x , UA + n y ^ A ) - (£+ fe)(n XiWB + n y , w J 

+ l -{ qi \Xx){XX\at !UA +q 2 \G)(Xx\4 !UJB + q 3 \Xy)(XX\al <UA + g 4 |G)(Xy| a t w t _ H .c), (5) 



where 2£ = u>\ — ll>2 is the biexciton shift, A = u>s — u>\ — 
LU2 — ^4 is the doublet splitting due to dot asymmetry, 
5a = <-^A — is the detuning between cavity mode A and 
transition frequency u)\ and 8 b = ^2 — lob is the detun- 
ing between transition frequency o>2 and cavity m ode B . 



These frequencies are shown schematically in Fig. 1(b) 



We now define a "balanced cavity" to be one for which 
the two cavity modes fall directly in between each of the 
doublets (8 a = 8b = A/2) and the exciton-cavity cou- 
pling constants are matched {q\ = q% and qi — q^j. An 
"unbalanced cavity" is one for which 8a,b 7^ A/2, and 
"unbalanced coupling" means that q\ ^ q% or q2 ^ q/±). 
We will show later that a balanced cavity accomplishes 
the required "which-path" erasure. 

The dynamics of states under the action of the time 
evolution operator, e~ lHt , generated by the Hamiltonian 
H is closed in the 12-dimensional space spanned by the 
basis B, which is shown in Fig. |^. 



Finally, we assume that the initial state of the system 
is biexcitonic, |^(0)) = |XX)|00)|00). 

III. DERIVATION OF MASTER EQUATION 



The theory of open quantum systems has been well 
studied (see e.g. |16|, [17)), and we adopt this formalism 
to analyse the exciton-cavity system interacting with the 
external continuum modes and measurement devices out- 
side the cavity. 

Wiseman |]l6f gives an expression for the master equa- 
tion for the conditional density matrix, p c , for a single 
measurement channel using imperfect detectors, while 
Gardiner and Zoller |l7| give a similar expression for 
many channels with perfect detection on each channel. 
Generalizing these results to an n channel conditional 
master equation with arbitrary efficiency detectors on 
each channel results in the conditional master equation 



|G)|10)|10)A|Xx)|10)|00)A|XX)|00)|00)A|Xy)|01)|00)A|G)|01)|01) 
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FIG. 2: The basis B for the evolution of the system along with transitions generated by the Hamiltonian H and cavity leakage. 
Coupling strengths between states are indicated. The top line of states spans the subspace of two excitations (i.e. exciton 
number plus photon number), the middle line spans the subspace of one excitation, and the single state on the last line spans 
the subspace of zero excitations. 



dp c = -i [H, p c ] dt + ^ I (rjj Tr { J jPc \ p c + (1 - r)j)Jjp c - ^(c]cjp c + p c c]cj 



dt 



J ( 3PC : - Pc ) dNj 



(6) 



r 



where H is the interaction Hamiltonian for the system, Cj 
is the system operator through which the system couples 
to channel j, Jjp c = Cjp c Cj is the jump operator for 
channel j, rjj is the detection efficiency of jump processes 
on channel j, and dNj is the jump increment. For the 
case where rjj = 1 for all j, this equation reproduces 
the result in Gardiner and Zoller |l7], §11.3.8.d], and for 
n = 1 it reproduces the result of Wiseman §4.1.2]. 

For the biexciton decay, there are several baths with 
which the exciton-cavity system is coupled. Firstly, the 
cavity modes decay at a rate k in order to couple the 
photons generated in the emission process to the outside 
world. This decay mode is coupled to the four cavity 
modes, where the system coupling operators are c\ — 



>c 2 



C3 



= \fna^ UJB and C4 = ^/na 



y,uj B ■ 



To quantify the effect of the cavity in erasing the fre- 
quency information, in what will follow, these channels 
are assumed to be perfectly detected, rji = ■ ■ ■ = r/ 4 = 1. 

Secondly, there may be spontaneous emission into pho- 
ton modes apart from those of the cavity. This decay 
channel couples via similar system operators, but with 
different decay rates, so that C5 = \ZT^\Xx)(XX\,cq = 
Vr~ s \Xy)(XX\,c 7 = VTT|G)(Xx| and c 8 = VTl|G)'(Xy|. 
These channels are considered to be inaccessible to an 
observer, so we set the detection efficiency to zero, 774 = 
■ • • = 7?8 = 0- For later sections, we will refer to these 
channels as "leakage channels" . 

Finally, we will add a phenomenological dephasing 
acting on the two exciton states \Xx) and \Xy). This 
is to simulate the effect of some unspecified bath (e.g. 
phonons) that is able to distinguish the intermediate ex- 
citonic state during the decay process. The system op- 
erators to which this bath couples are assumed to be 
c 9 = VT\i\Xx){Xx\,Cio = y / T~ l \Xy)(Xy\. Again, these 
channels are inaccessible to observer, so the detection ef- 
ficiency is zero, 779 = ?/io = 0. 

Since rjj — for channels 5 through 10, from | |T6| ] we 
have E[dNj(t)] = rjj Tr {Jjpjf)} dt = for j = 5 . . . 10, 



and since dNj(t) is non-negative, dNj(t) = for j = 
5 . . . 10. In accordance with the assumptions regarding 
channel efficiencies made above, the conditional master- 
equation between photon detections (i.e. dNj = 0, j = 
1 ... 4) , becomes 

4 

Pc= ~ i [H, p c ] + ^2 Tr {Jj ; p c } p c 
3=1 

10 10 ^ 

+ ^2jjPc-^2»(ctcjp c + Pcc}cj). (7) 



3=5 



3 = 1 



The second term in this equation is non- linear in p Cl re- 
flecting the fact that the evolution is conditional on the 
system not emitting a photon. For computational pur- 
poses, we convert Eq. (m into an equivalent linear equa- 
tion for an unnormalized density matrix p by defining 
Pc{t) — f(t)p(t), where f(t) is a scalar function to be 
determined. Substituting this into Eq. (^) gives 



fp + fp= - if{H,p}+f 2 Y,^{J 3 P}P 
10 10 



3=5 



3=1 



Collecting terms that are proportional to / and requiring 
that others vanish gives the linear, unnormalized semi- 
conditional (i.e. conditioned on only a subset, j = 1 . . . 4, 
of the channels) master equation 



p = -i [h, p] + J2 ftp - £ \i c h~p + p4<h), (9) 

3=5 3=1 



^y-j-jr ■ r^j^3 

along with the constraint equation for / 



fp=fJ2^{J3P}p- 

3 = 1 



(10) 
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FIG. 3: Schematic of a polarisation sensitive interferometer. 
A polarising beam splitter (PBS) splits the beam into x- and 
y-polarised paths and a relative phase cf> is added to one path. 
The paths are recombined, and the polarisations are rotated 
by 7r/4 using a half- wave plate (HWP), then detected with 
polarisation-sensitive single-photon detectors 



This can 

rt ,,x-^4 



be integrated to give / = 
{f : di£* =1 Tr {Jjp})' 1 Taking the trace of Eq. (g) 
gives Tr {p} = — Ylj=i Tr{J-/p}, an d so we see that 
/ = Tr{/5} , which is just the normalization condition 
for p c , i.e. p c = p/Tt{p] as required. 



intuitively clear from these two examples that the visi- 
bility of the interference fringes is an operational mea- 
sure of both the purity and entanglement of the input 
two-photon state, and is dependent on the off-diagonal 
elements of the density matrix, which are zero for non- 
entangled or completely mixed states. 

Interferometric methods for estimating ent ang lement 
have been discussed by Ekert and Horodecki p9| . They 
argue that d 2 — 1 separate types of interferometric exper- 
iment are required to estimate the entanglement of a pair 
of particles, d being the dimension of the Hilbert space 
for each particle. For the case of interest to us d = 2, 
so we expect three parameters will be sufficient to place 
bounds on the entanglement of the photon pair. In fact, 
we assume that anti-correlated states such as |xy) are 
never produced, which is a roughly consistent with ex- 
perimental observations showing that inter-exciton tran- 
sitions are rare, ]Tl[ , and so the number of experiments 
required is reduced to one. That is, we only need to 
measure a single visibility fringe in order to quantify the 
two-photon entanglement. 



IV. QUANTIFYING TWO-PHOTON 
ENTANGLEMENT 

We now develop a measure of the performance of the 
cavity in erasing "which-path" information. A polarisa- 
tion entangled photon pair is an archetypal example of 
a two-qubit system. Such bipartite systems have been 
studied extensively Q], and in particular, the entangle- 
ment of such pure bipartite system is well quantified by 
the von Neumann entropy of one subsystem. 

In the system we are concerned with, there is some sub- 
tlety however, since although the system-plus-continuum 
evolves to a two photon state, this is only determined 
once a measurement has been performed, observing both 
photons. Before the measurement, the system evolves in 
a much larger Hilbert space, so it is not entirely triv- 
ial to adapt measures like entropy to the case of interest 
here, not least because under certain circumstances the 
photon-pair is described by a mixed state. Instead we use 
an operational measure - the visibility. This arises natu- 
rally by considering the result of two-photon coincidence 
counting at the output of a polarisation sensitive interfer- 
ometer, depicted in Fig. ^, through which the photon-pair 
is directed. 

It is straightforward to show that a pure, entangled 
state of the form |xx) + |yy) passing through such an in- 
terferometer with a (f> phase-shift (per photon) on one 
arm will exhibit interference fringes in the two-photon 
coincidence counts between the detectors, with the co- 
incidence count rate proportional to 1 — cos(20). The 
factor of 2 in the argument of the cosine is a direct man- 
ifestation of the two particle nature of the state, and this 
has been observed experimentally |Q . 

Conversely, neither a completely mixed state such as 
|xx)(xx| + |yy)(yy| nor a pure, non-entangled state like 
|xx) will display interference fringes as <j> ls varied. It is 



A. Interferometry 

Quantitatively, we relate the output continuum field 
annihilation operators of the half wave plate, b' XLiJ ,b yjjJ , 
to the interferometer input field operators, & x ,w>&y,a>> ac- 
cording to 





i 






b x ,uj 




" 


-1 e -i * 







(11) 



We adopt the notation that a prime on an operator in- 
dicates that it is transformed consistently with Eq. pi]) , 
for example h X UJ = {b' x ^b' XUJ . 

The expectation of a two-photon coincidence measure- 
ment by detectors 1 and 2 will in general be given by the 
normally ordered (denoted by :...:) two-time correla- 
tion function (: n[ UB (t B )n' ]UA (t A ) :) c , where i,j € {x,y} 
itujj [I?], ^(J. The subscript c denotes the fact that the 
expectation is conditioned on the system having emitted 
zero photons in the interval [0, minji^, £_b}), and the or- 
dering of the operators in the correlation function will 
depend on the ordering of t& and t b ■ 

We may relate the cavity field output operators for the 
continuum mode I to the cavity input operators and the 
internal cavity operator according to b out (l, t) — b ln (l, t) — 
y/nai (t) . We will also assume that the cavity input is the 
vacuum so (bl n (l,t)b in (l,t')) = |fl7[] . Thus normally or- 
dered expectations of continuum modes may be replaced 
by normal- and-time ordered expectations of internal cav- 
ity modes, multiplied by a suitable power of ^/k. More 
detailed discussion of this point is given in Gardiner & 
Zoller |TJ. 

For example, the conditional expectation of detect- 
ing consecutive photons at detector 1 will be given by 
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( : "L^bI^aC^ and if *B > then B. Visibility 



Since the transformed operators in Eq. ( |T^ ) depend 



= 1 , r{i7x jta , s (*B)T(i B ,i > i){i7 x ' fc , A (<A)/'c(*>i)}} i (12) on according to Eq. ([TT]), we see that the quantity 

-K 2 Tr((n' tfa' TdviAtn' n(tA(n' )H\ ( : ™*-"b ^b)K^ a (*) : > must also depend on 0. Many of 

- « ^ u<W a x,. B j ^^Ajiax^Pc^jK-J li . the cro B ss terms ; anisll5 leaving the result 

where T(t B ,t A ) is the time evolution operator, which 

evolves the system from time t A to time t B , and for open (: n' x U)B (t B )n' x U)A (t A ) :) = k 2 (x + y + e 2l ^z + e~ 2l ^z*), 
systems it is non-unitary |jl7| . Very similar expressions (13) 

may be derived for the case where t B < t A . where 

I 



x = ( : a l,u B a^ B a[^ A a^ A :) = Tr {at >U)B a X;U , B T{a XjU>A p c at tU)A }} E [0,1] 
V = ('■ a l,u B a y ,uj B al UJA ay^ A :) = Tr {a\^ B a )/ ^ B T{a y . M)A p c al^ A }} € [0,1] 

z = ( : <4, UB a*,VB°y, UA °>x,<»A : > = ^{ a l^B a ^B T {a^ A p c al^ A }} € C, (14) 



which all depend on t A and t B though this notation has 
been dropped for brevity. We can define the visibility, V, 
from this expression to be the amplitude of the interfer- 
ence fringes divided by the mean (averaged over <j)) and 
it is 



V(t A ,t B ) = 



2M _ 
x + y' 



(15) 



Conceptually, V is the visibility of fringes generated 
by post-selecting photon pairs that arrive within the 
two-time window (t A ,t A + dt) x (t B ,t B + dt) as 4> 
varies. We note that we may compute the visibil- 
ity directly from p by making the definition x 



Tr{c 



X,U A PO'x,W A 



}}, with similar definitions 



for y and z, so that an equivalent expression for the vis- 
ibility is 



V(t A ,t B ) 



2151 



(16) 



C. Probability density 



We may also compute the joint probability density, V, 
for detecting a photon pair within the two-time window 



(t A ,t A + dt) x (t B ,t B + dt), as given in 0, §11.3.7 (d)], 



V(t A ,t B ) = 




x 1- 



Ia 



drJ2^nr,0)p(0) 
= E< : "U (W L u A (*a) :> c Tr {p{t A )} , 

= K 2 j2^-{x^ B (tB)T(t B ,t A ){j^ A (t A )p(t A )}} , 
ij 

= K 2 (x + y) (17) 



where we have again assumed t A < t B , although simi- 
lar expressions may easily be derived for t A > t B . The 
first factor in in the first equality is just the conditional 
probability density for either detector to register at times 
t A and t B given no emission beforehand, and the second 
factor is the probability of emitting zero photons in the 
interval [0, t A ). The second equality follows from Eq. ( |l0| ) 
and later. The third equality follows from Eq. (12) and 
recalling the fact that p c (t) — p(t) / Tr {p(t)} . Finally, 
Eq. ( |l7j ) shows that the probability density does not de- 
pend on <j) - detecting a photon pair after the interferom- 
eter occurs with the same probability density as detecting 
a photon pair before the interferometer, as expected. 
We also define the quantity 



P : 



V(t A , t B )dt A dt B . 



(18) 



In the presence of spontaneous emission into non-cavity 
photon modes, P < 1, indicating that not all biexciton 
decay events will be detected by the photodetectors fol- 
lowing the interferometer. We therefore interpret P as 



7 



the reduction factor of the two-photon detection rate, as 
compared with the biexciton pumping rate. 



D. Mean visibility 

We now define the mean visibility, which is a figure of 
merit for the degree of entanglement between the photon 
pair, 



•OO /*OG 

V(t A ,t B )V{t A ,t B )dt A dt B , 
o Jo 

2 roc />oo 



2k f f 
-5- / / \z{t A ,t B )\dt A dt B , 



(19) 



where we have divided by P so as to only count those 
decay events that are detected through the interferom- 
eter. If the visibility is unity (i.e. perfect erasure of 
"which-path" information), then V = 1, since the proba- 
bility density is normalised by P. On the other hand, if 
V is less than unity, then V will be also, so performing a 
two-photon interference experiment with all photon pairs 
produced will result in fringes of visibility V < 1. 

We see from Eqs. @, @ and @ that the quanti- 
ties we are interested in may all be determined directly 
from p, which makes calculations we perform in following 
sections simpler. 



E. Phase accumulation 

Whilst the visibility is a very important measure of 
the success of the scheme, since the initial state of the 
system, |Xx)|00)|00), is not an energy eigenstate, during 
the emission process, phase will accumulate at different 
rates on the xx- and yy-decay paths. The phase differ- 
ence accumulated between each decay path depends on 
the emission times of the two photons and is given by 
if> = axg{z(t A , t B )} = axg{z(t A , t B )}, corresponding to 
emission of a state of the form |xx) + e J¥ ' < '* A '* B '|yy). Since, 
for a given apparatus, ip depends only on the emission 
times t A and t B , this may be calibrated or computed, 
and hence accounted for, before an interference experi- 
ment (or whatever else is intended for the output photon 
pair) is done. If this phase is ignored, then the mean 
visibility will be lower than V since the description of the 
phase-averaged state will be mixed. 



F. Relation to Bell-inequality violations 

Instead of passing the photon pair through an inter- 
ferometer, we could imagine using an ensemble of such 
states to measure Bell- inequality violations. In particu- 
lar, we consider violations of a Clauser-Horne-Shimony- 
Holt (CHSH) inequality, where each photon is measured 
in one of two non-orthogonal bases specified by the angles 
9a and 9' A for the photon at of frequency uj a and B B and 




FIG. 4: Schematic setup for performing a CHSH Bell- 
inequality measurement. Photons of frequency lo a and lob are 
split with a dichroic mirror (DCM) then travel along paths A 
and B respectively, and are measured by rotated polarisation- 
sensitive detectors. 



9' B for the photon at of frequency u B [|Tj, as depicted in 
Fig. g. 

In terms of mode operators, the CHSH inequality re- 
quires the knowledge of correlation functions of the form 



E{6 



A,VB 



(d\d+ - d f _d-){c f +C+ 



c f c_) 



(: (dt d+ + dld-){clc+ + c f _c_) :) ' 



(20) 



where photon mode annihilation operators c and d are 
defined as 

c + = sin(0A)%,w A +cos(6 A )a x<UA , 
C- = COs(9 A )a y>ulA -sin^^Ox,^, 
d + = sm(9 B )a y ^ B + cos(9 B )a x ^ B , 
d- = cos(8 B )a yiUJB - sm(9 B )a XtUJB . 

We have not explicitly included times in these expres- 
sions, but we note that operators a\^ A act at time tA 
and a- uu/B at t B . It is straightforward to show that 



C + C + + C'_C- — Qy )W A Cty,u A + a x,U> A ^i^A ) 



~ a y,LJ A a Y,^A + a x,l 

c + c +-c- c - = cos ( 2 ^)(aJ,^ a y^A +ajU A °W) 
+ sin(2^)( 

with similar results for d±. As mentioned earlier, 
many cross-terms in the numerator and denominator of 
Eq. ( |20| ) cancel for the physical situation we consider, e.g. 
( : a t,u A a ^ u A a y,uB a y^B '■) = 0; so we can write it as 



E(0a, b ) 

= cos(26U)cos(26> B ) + 



x + y 



sm(29 A )sm(29 B ), 



= cos(26U) cos(26» B ) + Vcos(^) sin(26U) sin(20 s )(21) 

where x, y and z are defined in Eq. (|l4|). 
We now define the quantity 

B = E(9 A , 9 B )-E(9 A , 9' B )+E(6' A , 9' B )+E(9' A ,9 B ), (22) 

which, for classically correlated states satisfies B < 2 
This inequality is violated by certain entangled states 
such as Bell states, which are emitted from the biexciton 
system. We note that B depends on t A and t B , since V 



8 



B 



8 
6 

4 - 








2 






V cos(y) 


8 / 

6 

4 




0.6 


0.8 1 



FIG. 5: The quantity B versus V. The lower curve is B for 
■d — 7r/8. The upper curve is the maximum value of B for the 
corresponding value of V, allowing ■& to change with V. 



inequality. At <d — 7r/8 Eq. ( p4| ) reduces to 

B = \/2(l + Vcos(p)) (25) 

which is plotted in Fig. || (lower curve). We also show the 
maximum value of B for each value of V (upper curve), 
allowing $ to vary. The upper curve crosses B = 2 
at Vcos(i^) w 0.316, whilst the lower curve crosses at 
Vcos{ip) =\/2-l« 0.414. 

From these results we see that V and i3 are very closely 
related quantities. Since B may be computed from V for 
arbitrary angles, we will base our computations on V, 
from which a reasonable estimate for the maximum value 
B may be evaluated using Eq. (p^). Finally, we define the 
quantity B as 



and ip does, but we leave out the explicit notation. From 
Eqs. ( pl| ) and (^2|) we derive a linear relationship between 
B and V given by 

B = cos(26U) (cos(26» B ) - cos(26^)) 

+ cos(26' A ) (cos(26» B ) + cos(26^)) 
+ Vcos(<^){ sin(2(9 A ) (sin(26>s) - sin(20^)) 
+ sin(26^) (sm(2d B ) + sin(20' B )) } (23) 

We consider a special choice of angles that maximally 
violate the CHSH inequality, 6 A = 0,Q B = &,6' A = 2d 
and 9' B — 3i9, and using Eq. ( pl| ) in the expression for B 
we find that 



B = cos(2i?) (3-2 cos(4tf ) + cos(8tf)) 
+ 8cos(2tf) 3 sin(2tf) 2 Vcos(^). 



(24) 



For V = 1, this gives B = 3cos(2z9) — cos(6z9), which has 
a maximum at •& = vr/8 of 2V2 > 2, violating the CHSH 



B 



1 

P 



B(t A ,t B )T(t A ,t B )dt A dt B , (26) 



where we have again divided by P in order to count only 
those photon pairs that are detected in the experiment. 
G. Phase-averaged Bell-inequality violation 



Comparing B with Eq. ( B4j ) or Eq. (25), we see that 
we need to compute the quantity 



1 

P 



W cos(tp)dt A dt B . 



We generalise this to account for the possibility of adding 
a fixed relative phase <f> to one decay path (e.g by adding 
a phase plate on the y-polarised photon path, as in the 
interferometer stage of Fig. ||), so that cos(<p) — ► cos(</? + 
(f)). We maximise the above integral over <p to arrive at 
the phase- averaged visibility 




W cos(ip)dt A dti 



VV s'm(ip)dt A dti 



(27) 



The phase-averaged visibility, Q, gives the visibility 
of the fringes in a two-photon interference experiment 
where no attempt is made to resolve the phase accumu- 
lation. In regard to a CHSH-inequality violation exper- 
iment without sufficiently fast time-resolved detection, 
violations may still be seen if Q > 0.316, since the func- 
tional relationship between Q and B is the same as that 
between B and Vcos(<p), as shown in Fig. |H[ 



V. ANALYTICAL RESULT FOR BALANCED 
CAVITY 



We now show that for a balanced cavity (i.e. S A = 
5b = A/2, qi = 53 and (72 = 94), in the absence of spon- 
taneous emission and dephasingjE^ = T s = 0, the model 
predicts that the visibility is unity for all (t A ,t B ). Since 
we assume = T s = 0, we may write the unnormalised 
density matrix as p(t) = \^(t)){ip(t)\ and then Eq. (|^) 
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may be written as a Schrodinger equation 

j i \m) = -m\m) 



(28) 



for the state vector, \ip(t)) , with a non-Hermitian effective 
Hamiltonian given by 



H = H - in/2(h Xil 



+ hy !UB ). (29) 



We do this by considering the transformation of the ef- 
fective Hamiltonian and state vector under exchange of 
the polarisation, \Xx) <-> \Xy) and |01) <-> 1 10) . This 
transformation, denoted hereafter by is just a permu- 
tation on the basis elements, leaving the two elements 
|XX)|00)|00) and |G)|00)|00) invariant. A matrix repre- 
sentation of # shows that is both orthogonal and sym- 
metric. 

Consider evolution in S 2 - Swapping x and y polarisa- 

# 



We write the solution to Eq. (§|) as = e~ mt |0(O)). tions ma P s U 2 ^ H 2 = -H* 2 and |0(t)) 2 \xp(t)) 

As a result the time evolution operator (e l " 2 *)# 



The smooth evolution is of course punctuated by quan- 
tum jumps, corresponding to photon detections following 
the interferometer. 

Reformulating the equations of motion in terms of 
quantum trajectories Jf?} has several advantages, and 
most obviously it reduces the number of unknown quan- 
tities, since we can now solve for the state vector rather 
than the density matrix. It is straightforward to show 
that the effective Hamiltonian, Ti, only couples states 
within the same excitation-number subspace. Coupling 
between the zero, one and two excitation subspaces (de- 
noted So, Si and 52 respectively) occurs only during the 
jumps, and the excitation number irreversibly decreases 
by one at each jump as photons leave the cavity. Thus, 
for the smooth evolution between jumps, we may con- 
sider the evolution restricted to states within each Sj 
independently, and for each Sj we consider the effective 
Hamiltonian, Tij, restricted to that subspace and acting 
on the state vector \ip(t))j. 

Using the quantum trajectories formalism, we find 



X = (^(*.A,*b)|$x(*A>*b)) 
V = Wv(*A,*B)|</>y(*A,*B)) 
Z = {$ y (tA,tB)\'<Px{tA,tB)), 

where, assuming t A < t B = t A + T we have defined 



(30a) 
(30b) 
(30c) 



\i>\{tA,t B )) = a\^ B e 



— %H\T 



MtA)) 



= a- liUB e- iH ^a- hUA e- iH ^\i,(0)}. (31) 

A very similar expression exists for tA > t B , and the 
following reasoning applies equally to both cases. The 
state vector \ip\{tA,t B )) £ So, since the initial condition 
1-0(0)) = |AA)|00)|00) e S 2 and the effect of the two an- 
nihilation operators in Eq. ( |3l| ) is to reduce the excitation 
number by two. 

The one dimensional subspace So is spanned by the 
system ground state |G)|00)|00), so a state |V>(f)) £ So 1S 
mapped smoothly to a scalar, ip(t), by the trivial map- 
ping ee ((G|(OO|(OO|)|0(i)). We may therefore write 
x, y and z in terms of the scalar quantities ip x (t A , t B ) and 
Tp y (t A , t B ): x = 0*0 X , V = rpy rpy and z = V>*0 X , where we 
have dropped the time dependent notation for clarity. 

In what follows, we establish that for a balanced cavity, 
X and y are related by a unitary factor. This means 
that they have the same amplitude, from which it fol- 
lows that the visibility is unity for a balanced cavity. 



e in * t — e in2t when acting on states in 52- We note 
in passing that for an unbalanced cavity/coupling H.^ ^= 
—Ti-2, which is why it is critical that the cavity be bal- 
anced this argument to be valid. 

A similar result applies to evolution in Si , except that 
the effective Hamiltonian "Hi does not transform under 
polarisation swapping quite as simply. Instead, it may 
be shown that Hi - H d -► 7if - H* = -(HI - H d ), 
where H d is a Hermitian matrix acting on elements of Si 
and satisfies [Hi,H d ] = 0. Thus Ti* = -H\ + H d + H* 
and {e~ m ^)* = e" iW ?* = ^nit e -r{H d +H#)t _ The fac _ 
tor U d (i) = e~ 4 (^+ff*)* j s un itary, since H d is Her- 
mitian. In particular, U d (t) acts on states of the form 
o.\^ A e~ mtA 1^(0)) € Si in a simple way: it multiples the 
state by a time dependent unitary scalar, e l6t . 

Having established the effect of # on the time evolution 
operator acting on Si and 52 we see that for example 



\^{tA,t B ))* = (a x , UB e- iHT a : 



a e- m * T a e - m * tA 



l^(0))) # , 



1^(0)), 



iH'r, 



,-tW-r, 



jn*t A 



e WT {a. 
e i6r $ y (tA,tB)Y 



10(0)), 
1^(0)))*, 
(32) 



This second 
|XA)|00)|00) # 



line 



follows since 
|AA)|00)|00) = 



|0(O))# = 
|0(O)) and 

a x ,oj# = a y ,uj and the third line follows by consid- 
ering the arguments in the preceding two paragraphs. 

On the other hand, since \ip x (tA, t B )) £ So it is evi- 
dent that \^ x {t A ,t B ))# = \i> x {t A ,t B )) as |G>|00)|00) is 
invariant under Together with Eq. (]32]), this implies 
|0x(*A)*s)) = e l6 ' T |0y(i J 4, t B ))*, and we conclude that 
0x = e l£,T 0*. It follows that x — 0*0 X = V>y0y = y, and 

also \z\ — |0*0U = |e~ l6,r 0xl = x. Using these two re- 
sults and Eq. (16) we see immediately that for a balanced 
cavity V = 1, proving that the visibility is unity for all 
times. 



VI. NUMERICAL RESULTS 

In this section, we present the results of computations 
for V for unbalanced systems, and results for Q which 



10 




FIG. 6: In each row, from left to right, k = 0.5,1,2,4 and Td — T 3 — in all figures, (a) Time dependent probability 
distribution, V(tA,tB), for a balanced system, (b) Visibility, V(tA,ts), for an unbalanced system, where 91,2,3 = 1 and 
Q4 = 1.1. (c) Relative phase, <p(tA,tB), for a balanced system. Superimposed on each panel is a contour plot of V '. 



characterizes the extent to which Bell-inequality viola- 
tions may be observed. For the problem parameters, we 
take experimentally relevant values typical for GaAs self- 
assembled dots A = 50 /ieV @, f = 0.5 meV, q = 50 
fxeV |2^, K > 100^eV, though values of k much lower 
than this may be possible with novel hemsipherical cavi- 
ties pfj. Throughout this section we rescale all energies 
so that qi = 1, A = 1, £ = 10. Time is also rescaled 
accordingly, so one time unit corresponds to 83 ps. 



Figure 6(a) shows plots of the probability distribution 
of emission times for a balanced cavity with no leak- 
age channels. Numerically computed visibility is unity 
to within numerical accuracy and V = 1 to within 
10~ 4 , when integrating out to tA = ts = 200. Notice 
Rabi oscillations in emission time for strong coupling 
(/t < qi = 1) and exponential decay for weak coupling 
(« > q). For strong coupling there is a significant prob- 
ability of emitting photons in either order, but in weak 
coupling the order ts > tA is strongly favoured, indicated 
by the sharp edge along tA = is- 



We also note that in the weak coupling regime, V has 
a tendency to broaden with increasing k, which is some- 
what counter-intuitive, since larger n corresponds to a 
more leaky cavity, and one would expect the photon- 
component of the internal state to leak away more 
rapidly. However, this phenomenon may also be seen 
in the much simpler case of a single two-level atom inter- 
acting with coupling rate q with a single optical mode of 
a leaky cavity. In that case it is straightforward to show 
that there is in eigenvalue of the effective Hamiltonian 
for the open system given by q 2 KT x /2 + 0(k~ 2 ), which 
corresponds to a long time constant for large k. When 
k w q there is a kind of impedance matching, and the 
temporal extent of V is smallest. 

As established previously, the visibility is unity for a 
balanced system. For an unbalanced syste m the visibility 
drops below unity, as shown in Fig. 6(b) where (74 = 1.1 
(with 91,2,3 = 1), for different values of k. The probabil- 
ity density V for this case looks very similar to 5(a)| so 
is not shown here. We note that the visibility depends 
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(c) (d) 



FIG. 7: Dependence of (a) V and (b) Q on 34 for various n. Dependence of (c) V and (d) Q on 8 a and 8b- V and Q are not 
sensitive to the sign of 8a or 8b so other quadrants look similar and are not displayed. The grey plane at Q — 0.316 demarks 
the threshold above which Bell-inequality violations may be observed. 



only on ts when tA > tg, i.e. it is frozen at the value it 
reaches at tg- Notice that the probability, V , of emitting 
a photon pair is small at the the same times that V has 
large excursions from unity, which means V is not affected 
as much as one might expect given the large fluctuations 
in V. The difference between strong and weak coupling is 
striking, again with oscillations being replaced by decay. 



which shows phase singularities. 



The relative phase, <p, is shown in Fig. |6(c)| for a bal- 
anced system. Also superimposed on each panel is a con- 
tour plot of the emission probability density, V . In the 
strong coupling regime, during Rabi oscillation peaks, the 
phase accumulates relatively slowly, with rapid phase ro- 
tations in between. In the weak coupling regime, the 
phase accumulates at a fairly constant rate, which is 
roughly proportional to A. The diagonal stripes indicates 
that in weakly coupled cavities the phase accumulation 
depends only on the time interval between photon emis- 
sion, tg — tA, in contrast to the much more complicated 
dependence of the phase in the strong coupling regime, 



Figure 7(a) shows the variation of V versus (74 for var- 
ious k. For = one decay path is turned off so we 
expect completely non-entangled photon pairs, and this 
is evident in Fig. 7(a) as V = when q^ = 0. We also 
expect that V = 1 when 54 = 1, since then the couplings 
are again balanced, and this also is evident in Fig. [7(a) . 
The variation with qi is identical to that displayed here, 
whilst the variation of V with q\ and qs is qualitatively 
very similar, so is not shown here. 

As discussed earlier, Q is a significant quantity which 
determines whether the photon pair can produce Bell- 
inequality violations in the absence of time-resolved de- 
tection, so that phase is ignored. In particular, as shown 
in Fig. H if Q > 0.316 then the photon pair can produce 
Bell-inequality violations, even in the case that ip is ig- 
nored. Figure 7(b) shows Q for the same values of 174 
and k as in Fig. 7(a), where the grey plane demarks the 
threshold, Q = 0.316, to see Bell-inequality violation, as 
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FIG. 8: Q versus A for k. Varying A changes the rate 
of phase accumulation between the photon detection events. 
Grey plane as in Fig. 



it will in all following plots of Q. Whilst Q is everywhere 
less than unity, there are parameter values where Bell- 
inequality violation may still be observed without using 
time-resolved detection. 

If the cavity geometry is such that <5j4^A/2or#B^ 
A/2, then the cavity is unbalanced. This has the effect of 

Results are 



reducing V below unity, as shown in Fig. 7(c 



only displayed for Sa.b > A/2, but the other quadrants 
are similar. Figure 7(d) shows Q for the same parameter 
values, again with the plane denoting the threshold for 
Bell-inequality violation. 

The phase accumulates in between the photon detec- 
tion events roughly at a rate proportional to A, the split- 
ting due to dot asymmetry, and this directly affects Q, 
since for smaller A we expect the phase to be more nearly 
constant over the photon emission lifetime. This may be 
seen in Fig. |8|, where for small A, Q approaches unity, 
although Bell-inequality violation may still be seen for a 
wide range of k and A. 

The two different leakage channels that we consider 
in this paper are spontaneous emission into non-cavity 
modes, which occurs at a rate r s , and dephasing which 
happens at a rate Td- 

The spontaneous emission does not affect the visibility 
of those photon pairs that arrive at the detector, but it 
does change the rate of detection, since some photons are 
lost. The reducti on in photon detection rate, given by P, 
is shown in Fig. |9(a) . The roll-off in P is roughly pro- 
1 Surprisingly, the spont aneou s emission 

although 



portional to T s 
enhances Q, which may be seen in Fig. 9(b) 



we note that the source is then no longer deterministic. 
Experimental work indicates that high /3-factors, up to 
f3 = 0.9, are possible and recently (3 — 0.83 has been 
observed p5[ , which is the fraction of photons emitted 
into the cavity mode, corresponding to the fraction of 
photons emitted into the desired cavity modes. Inter- 
preting P as the /3-factor, from Fig. 9(a) we surmise that 
for P 0.9, the experimentally relevant range of the 



spontaneous emission rate is T s <C 0.1 = 5 /ieV, which is 
a regime in which spontaneous emission is negligible. 

Figures 9(c) and |9(d)| shows the effect of the phe- 
nomenological dephasing term, Td, for different values 
of k. V and Q decay roughly as T^" 1 . In all panels of 
Fig. ^, there is a peak along n « q which is due to the 
fact that V is temporally narrowest when this condition 
is met, and hence there is less time for leakage to take 
place. For very cold temperatures, around 1 K or lower, 
pure dephasing rates have been observed to be around 
1/ieV ^6j, corresponding to = 0.02, which is negli- 
gible. For higher temperatures, the pure dephasing has 
been observed to in crease at roughly 0.5 - 1.6 ^teV/K 
H, [27j. From Fig. gd|, the pure dephasing becomes 
important near < 1 corresponding to a temperature 
between 30 and 100 K for the experimentally relevant 
range given above. 



VII. DISCUSSION 

In the previous section we found that the numerical 
results for a balanced system concur with the analytic 
result derived in Section [v| where we established that 
the visibility is unity in this case. We also noted that 
V is degraded by any effect which may cause the cav- 
ity or couplings to be unbalanced. Imperfections in the 
cavity geometry will result in an unbalanced cavity so 
5a,b 7^ A/2, and it was shown above that this reduces 
V. Similarly unbalanced coupling constants also results 
in decreased visibility. 

Both of these effects may be understood heuristically 
using a much simpl er mo del whic h captures the gross 

Firstly, we note that 



features seen in Figs. 6(b) and 7( 



a two-photon state given by a x |xx) + %|yy) will produce 
two-photon interference fringes with visibility 



V 



2\a T a v 



(33) 



Secondly, we make two ad hoc simplification s of t he level 
structure of the quantum-dot shown in Fig. 1(a). These 
simplifications are (i) to ignore the crystal ground state, 
|G), and the corresponding transitions thereto, and (ii) 
to treat the remaining three level system, composed of 
| XX), \Xx) and \Xy) as a pair of independent two-level 
systems (TLS), {\g)i, |e)i} and {\g)2, |e)2}, each of which 
interact with one of a pair of degenerate cavity modes 
distinguished by polarisation. With these two assump- 
tions, the energy level structure becomes that pictured 



in Fig. 10(a) 



The physical motivation for these seemingly arbitrary 
assumptions is firstly that once the biexcito n dec ay pro- 
ceeds along the x- or y-polarised paths of Fig. 1(a) , the re- 
sulting two-photon amplitudes, a XiY are determined, even 
though the dynamics of the emission are not complete. 
Thus the two-photon amplitudes are largely determined 
by the initial single-photon decay process, justifying (i). 
Secondly, whilst the sum of probabilities to take the x- 
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or y-polarisation decay paths is unity, apart from this 
constraint the rate equations for the two decay processes 
are otherwise uncoupled, so the system similar to a pair 
of uncoupled TLS', one for each decay path, justifying 
(ii). Ultimately, this highly simplified model will be ver- 
ified by its qualitative agreement with the more realistic 
model discussed throughout this paper, and its value is 
in the intuition it lends about the origin of the effects 
seen in the numerical calculations. 

A TLS interacting with a cavity mode is well under- 
stood in terms of the Jaynes-Cummings (JC) model [ fl4| . 
For TLS' initially in the state |e>< (i = 1, 2 is the TLS la- 
bel) with energy spacings Vi, oscillator frequency u>i and 
detuning 6j — V j — u>i, with TLS-cavity coupling rate f2j, 
(see Fig. 10(a)| ), the time-averaged photon population is 
given by pi = Q?/(2i2?), where Ri = y/6f + Hf, and we 
conclude that the average amplitude of photon occupa- 
tion satisfies 



*x,y| 



V2Ri. 



(34) 



We now compare the predictions of this simple model 
with the more complete one for an unbalanced cavity, 
wherein the cavity mode is not tuned to the mean of the 
transition frequencies, u> ^ (v± + v-i)j1. If the TLS' are 
detuned by an amount 6\,i = 6 =F D/2 respectively from 
the degenerate cavity modes, (see Fig. 10(a)), each with 
the same cavity-coupling strength Oi : 2 = ^, the visibility 
is then given by 



V 



2R1R2 



Rf 



R\ ' 



(35) 



where we have taken | a x . v | fr om Eq. 



). This expres- 
sion is plotted in Fig. 10(b) as a function of 6 (using 
il = 0.61,13 = 1) along with V (using (71,2,3,4 = 1,K = 
0.4, A = 1,6 A = A/2,C = 10,r d)S = 0). ' Clearly the 
forms of the two traces are in qualitative agreement 
demonstrating the heuristic validity of the simple model. 
The value ft = 0.61 is_ selected to fit Eq. @ to the 
numerically computed V, but it is of the same order as 
Qi = 1. 
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We also compare the predictions of the simple model 
for unbalanced coupling to the realistic model, and so we 
take fii 7^ 2 , but suppose the detunings between the 
two-level systems and their respective harmonic oscilla- 
tors are equal, Si = —62 = 5. It is straightforward to 
show that 



V = 



(i?i^i) 2 + (R 2 n 2 y 



(36) 



which is plotted in Fig. 10(c) (using = 0.42, 5 = 1.69 



both fitted para meters ), along with V (other parameter 
values as in Fig. 10(b) ). 

The simple, heuristic model of two uncoupled two-level 
systems predicts a visibility, V, that is qualitatively in 
agreement with V calculated using the complete model 
discussed in earlier sections. Thus we can understand the 
most significant effect of variation of qi and Si on V is to 
change the relative amplitudes to ta ke each of the two 
decay paths illustrated in Fig. 1(a). Since the photon 



pair is only maximally entangled when the amplitudes 
of the |xx) and |yy) components are equal in magnitude 
(i.e. for the state (|xx) + e~ ll ^\yy})/^/2), parameter varia- 
tions that result in unequal decay path amplitudes result 
in sub-maximally entangled photon pairs. Such parame- 
ter variations correspond directly to the situation of an 
unbalanced system. 

The analysis above gives us some further insight into 
the decay process. The maximum amplitude of the pho- 
ton excitation is f2 2 /i? 2 and so the leakage rate of photons 
from the cavity will be suppressed by this factor. That 
is, we expect that the rate of decay of excitation from 
the cavity will be roughly refi 2 /(fi 2 + 5 2 ). Therefore, as 
the detuning S increases, the photon emission rate slows 
roughly as ~ 1/S 2 for 6 > Q. This will mean that for 
detunings significantly larger than the coupling strength, 
leakage effects will become significant - the lifetime of 
the excitation in the cavity will become comparable to 
the decay rate for dephasing or spontaneous emission. 

It is also worth noting that when S ~ 2£ (i.e. the 
exciton-cavity detuning is near the biexciton shift) the 
model developed in section || breaks down, since signif- 
icant cross coupling between exciton states and cavity 
modes will set in. 

Spontaneous emission decreases the detection proba- 
bility, which could be corrected with post-selection, since 
only events in which two-photons are registered count 
towards the measurement, and as mentioned previously, 
experimental work has shown that this is negligible for 
experimentally relevant systems [^]. Dephasing of in- 
termediate states decreases the visibility exponentially in 
time. Temperatures of a few Kelvin provide sufficiently 
low dephasing rates that it is negligible, but the exci- 
tonic dephasing becomes important at temperatures of 
several tens of Kelvin, [^6[ [27|]. In principle, these effects 
may be distinguished using sufficiently fast time-resolved 
spectroscopy, since spontaneous emission will result in 
fewer photons reaching the detector, whereas dephasing 
would result in a time dependent visibility that degrades 



exponentially with time. 

So far we have not addressed the issue of how to ex- 
perimentally construct a cavity with the required spec- 
trum, shown in Fig 



1(a) 



and a detailed proposal for 
its implementation is beyond the scope of this paper. 
The enhanced exciton emission into the cavity mode is 
known as the Purcell effect and requires small cavity vol- 
umes, so that the exciton-cavity mode coupling strength 
is large and the density of available photon modes is small 
1^4] , |2§f| . Thus small cavities are necessary, and high Pur- 
cell factors have been demonstrated experimentally in 
single- wavelength sized cavities |Q ^9| . 

In contrast to the need for small cavities is the rela- 
tively small biexciton shift, 2£, which is around 1 meV. In 
order for a single Fabry-Perot resonator to accommodate 
modes spaced by 1 meV (i.e. the free spectral range, 
FSR), the cavity length would need to be of the order 
of 100/xm or more. This is too long for several reasons, 
because the Purcell factor decreases with cavity length, 
and also because growth of such a large heterostructure 
would be prohibitively difficult. As a result the cavity 
to which we have been referring throughout this paper 
would need be based on a more complicated geometry 
than merely a pair of distributed Bragg reflectors (DBR) 
forming a linear resonator. 

We stress that a more complex geometry is not just a 
requirement of this proposal, but that it would be neces- 
sary for a system even with symmetric quantum dots. If 
the cavity did not have separate modes near the exciton 
and biexciton doublet frequencies then only one tran- 
sition could couple strongly to the cavity, and the other 
transition would be sufficiently off resonance (5 > q) that 
the Purcell effect for that frequency would be suppressed, 
i.e. either the biexciton-exciton or exciton-ground tran- 
sitions may be well coupled to the cavity, but not both. 

It may be possible to engineer a small cavity with a pair 
of closely spaced modes using photonic crystals. If, dur- 
ing the growth of each DBR stack one layer was permit- 
ted to grow to be larger than A/4, then the cavity would 
look more like two coupled cavities, which may have the 
desired split modes. Certainly, geometric effects in mi- 
cropillars have been shown to produce a pair of modes 
spaced by ~ 5 meV p9| , though this was due to lift- 
ing polarisation degeneracy with elliptical cross-section 
cavities, which is undesirable for our scheme. 

Experiments using confocal hemispherical cavities, 
consisting of a planar Bragg reflector at the focal plane 
of a hemispherical reflector, of length 50 to 1000 jiva. are 
currently underway for quantum information processing 
purposes |3(| . In this configuration the cavity mode waist 
diameter is of comparable size to the optical wavelength 
and coincident with a quantum dot so that the exciton- 
cavity mode coupling strength is reasonably large. This 
arrangement may provide the two requirements of the 
present paper: both strong coupling between the dot ex- 
citations and the cavity mode (up to several tens of ^eV) 
and small FSR so that each doublet is on resonance with 
a nearby mode. It is quite plausible that by tuning the 
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1(b), thereby 



cavity length to vary the FSR and applying an external 
DC electric field to induce a Stark shift in the doublet 
frequencies, one may bring both doublets close to cavity 
modes simultaneously, as depicted in Fig. 
realising the requirements of this proposal. 



VIII. SUMMARY 

We have shown analytically that using a cavity with 
a particular mode structure facilitates the production of 
polarisation entangled photon pairs from an asymmetric 
quantum dot, which otherwise produces photon pairs en- 
tangled in both polarisation and frequency. We demon- 
strated this by computing the visibility of two-photon in- 
terference fringes produced using photons generated from 
such a cavity-quantum dot structure, and related this to 
their potential to demonstrate Bell-inequality violations. 

We have quantified the effect of various errors in the 
cavity mode structure, showing that the visibility is 
not degraded badly by mistuned cavities or unbalanced 
dipole coupling strengths, and for experimentally ac- 
cessible regimes is above the threshold at which Bell- 
inequality violations may be detected. 

Of major significance to this scheme is the phase 



accumulated between single photon emission events. 
By defining the phase-averaged visibility we were able 
to compute the effect of ignoring this phase on Bell- 
inequality and two-photon visibility measurements. Such 
phase-ignorance arises when the available time-resolution 
of the photon detection apparatus is longer than the 
asymmetry splitting, A. 

We showed that ignoring the phase reduces the ef- 
fective entanglement, but there are still experimentally 
accessible regions of parameter space that exhibit Bell- 
inequality violation, even when the phase is ignored, and 
the two-photon states are thus potentially useful sources 
of entanglement. 
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